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We discuss zero-frequency transport properties of various 
spin- 1/2 chains. We show, that a careful analysis of Quan- 
tum Monte-Carlo (QMC) data on the imaginary axis allows to 
distinguish between intrinsic ballistic and diffusive transport. 
We determine the Drude weight, current-relaxation life-time 
and the mean-free path for integrable and a non-integrable 
quantum-spin chain. We discuss, in addition, some phe- 
nomenological relations between various transport-coefficients 
and thermal response functions. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



I. INTRODUCTION 

The role of the spin excitations on the transport prop- 
erties of quasi-one dimensional Mott-insulators has been 
the subject of extensive experimental rescapch in the last 
few years. A recent 17 NMR investiffiationlil of SraCuOa, 
extending an earlier 63 Cu NMR studyB, measured a q = 
spin-diffusion coefficient (equivalent to diffusive magneti- 
zation transport) several orders of magnitude larger than 
the value for conventional diffusive systems. Thermal 
transport measurements in Sr2Cu03 and SrCu02 indi- 
cate at the same time, quasi-ballistic transport with a 
mean-free path of several thousands! 3 ] of A. _ 

It is well known from structural consideratiqeap and 
from studies of the magnetic excitation spectrumPfl, that 
Sr2CuC>3 and SrCuC>2 can be accurately described by the 
XXZ chain 



jj(xxz) 



J XX 



{sf S i+l 



S\ $1+1) 



Evidence for ballistic (or quasi-ballistic) magnetization 
transport have been found in recent exact djagonaliza- 
tion studiesQ of H^ xxz ^ at high temperatureaJ'H. A con- 
nection between integrability, conservation laws and bal- 
listiC|-tra.nsport has been proposed by Zotos and cowork- 
ersaOEI. If the current-current correlation does not de- 
cay to zero for long times, i.e when part of the current 
operator is conserved, i.e. when a certain (non-zero) pro- 
jection of the (here magnetic) current operator commutes 
with the hamiltonian, the transport is ballistic even at fi- 
nite temperatures. This seems to be the case, in general, 
for Bethe-Ansatz solvable models like H^ xxz \ although a 
formal proof for this connection is still outstanding. 

At present it is unclear, whether there exist non- 
integrable models which do exhibit ballistic transport, 
none has been found so far. Real compounds like 



Sr2Cu03 and SrCuC>2 correspond to H^ xxz ) anyhow only 
in first approximation. It is therefore important to ex- 
amine whether general, non-integrable, quantum spin 
chains show ballistic or diffusive transport properties 
at finite temperatures, r^jhis question has been stud- 
ied by Rosch and AndreiEj within a short-time approx- 
imation (a memory-matrix— approach, extending an ear- 
lier analysis by Giamarchillil) for Luttinger-liquids with 
higher-order Umklapp-scattering. They found only expo- 
nentially small deviations from ballistic transport away 
from commensurability. An alternative route to diffusive 
transport, the coupling spin-phonon coupling has been 
studied by NarozhnyPl . 

In his seminal papert£l in 1960, W. Kohn proved that 
the existence of a the delta peak at zero frequency (the 
Drude peak) in the conductivity is the essential difference 
between ground states with localized and extended elec- 
tronic states. A simple extension of this idea to the spin 
transport can be used to distinguish, without other ex- 
plicit information about the excitation spectrum, a spin 
insulator, like spin-Peierls compounds, from a spiB-con- 
ductor like Sr 2 Cu03. Despite an on-going efforttlUlS de- 
voted to this problem, the fundamental difference be- 
tween models with ballistic and diffusive transport prop.- 
erties has shown up only recently in QMC-simulationsE3. 
The purpose of this paper is to explain in detail how this 
important issue can be tackled numerically by QMC. 

The organization of the paper is as follows. Section || 
contains the basic definitions and sets the notation that 



we will use along the paper. In Section III we deduce the 



connection between the spin current-current and density- 
density correlation functions emphasizing the role of 
boundary terms that occur in Matsubara formalism and 
in Section ^ we discuss how to exploit that connection to 
compute the conductivity in imaginary frequency using 
Quantum Monte Carlo cluster algorithms in an efficient 
way. In Section [v| we describe a procedure to extract 
transport coefficients (Drude weight and the diffusion co- 
efficient) from QMC-data in general ID interacting sys 



terns either integrable or non-integrable. In Section VI 
we apply this method to the XXZ chain and we obtain 
the Drude weight at finite temperatures. We then dis- 
cuss several phenomenology relatio ns b etween tr anspo rt 
and thermal coefficients in Section VII. Section VIII is 



devoted to the computation of diffusion constants, mean- 
free paths and life-times in a non-integrable spin chain. 
In Section IX we present our conclusions. 
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II. SPIN CONDUCTIVITY 

QMC-simulations yield in general correlations func- 
tions on the imaginary-time axis. We therefore consider 
the Kubo formula for spin conductivity in the Matsubara 
formalism. 

The spin conductivity in one-dimensional spin chains 
can be defined as the response of the current to a ho- 
mogeneous and time dependent twist $ — (pi in the 
quantization axis. 

)($) = £ [hcos^i) + jfsinifr) + J Z S?S? +1 ] 



where 



and 



Jl — 2i \ 1 '+ 1 ; M-lJ 



(1) 



(2) 



(3) 



Formally (|l|) is the hamiltonian of a XXZ chain in which 
the quantization axis of the local spin operators has been 
rotated by a site-dependent angle <fti along the z-axis. To 
obtain the expression of the Kubo formula for the spin 
conductivity we expand H( xxz \<&) in a Taylor series: 



where (K) — QT\ ki) is the expectation value of the ki- 
netic energy per site and A is the current-current corre- 
lation as a function of the Matsubara frequency, 



e iUnT (j z (q,r)f(-q,0)) dr . (9) 



The response to the time-integrated twist is then ob- 
tained from (|8|) and the dynamical conductivity takes 
the usual form 

a[q ^ n) = -(K)-k{ q ^ n ) ^ D(q, Un ) (1Q) 



Eq. (|1C|) leads via analytical continuation iu n — > uj + iS 
and 5 —> 0, using a(ui n ) = liniq^o c(g, u> n ), to the usual 
representation of the dynamical conductivity 

a(w) = nD(T)6(cj) + a reg (cj,T) , (11) 

where D(T) is the Drude wcight0 that can be computed 
via D(T) = lim Wri _ (0 lim g ^ D(q,uj n ): 



D(T) = -(K) - A(g->0,a; n ^0) 



(12) 



It can be provecti that this limit4s consistent with the 
definition introduced by W. Kohrfi3 at T — 0, 



D 



(13) 



= + 3" <t>i + — 



Ij2 



(4) 



and we obtain the total spin-density current by differen- 
tiating with respect to ipi. 



dH {xxz \<5>) 



= tiff = 3? + hti 



(5) 



The first term is the paramagnetic part of the current. 
If the z-component of the magnetization is conserved, 
it can also be deduced using the discretized continuity 
equation. 



J^f(*) + (if(*)-if-i(*)) = 



(6) 



where the second term is the discrete version of the di- 
vergence in one dimension. If we combine it with the 
equation of motion 



0_ 

dt 



sm = i[H,s? 



(7) 



we obtain the expression (|3|). The second term in Eq. (§), 
proportional to the magnetic flux is called the diamag- 
netic current. The expectation value of the total current 
is 



(j T (q,u n ) = -( (K) + A(q,u n )) fa , 



(8) 



where E is the ground state energy and the total exter- 
nal flux. Recently D has been extended by Zotos, Naef 
and PrevlosekM to finite T 



D(T) = LJ2 



Z V d 2 $ 



(14) 



It is important to note that the limits lim g ^o and 
lim^^o do not commute. When the limits are taken in 
the opposite order one obtains the conventional spin stiff- 
ness which represents the response to a static twist. 



Ps 



lim lim D(q,u) n ) 

q— >0 UJ n ^0 



(15) 



A non-zero value of the Drude weight implies that 
the total magnetic-current does not decay to zero when 
t — > 00 ( i.e. the transport is ballistic). The most simple 
example to illustrate this situation is the XX chain. In 
that case [Hxx, j z {q)} = 0. Taking the spectral repre- 
sentation of A(g, uj n ) for uj n > we have 



Mq,u n ) = — ^ 

m.7i 



1 e-^\(m\j*(q)\n)\ 2 



IU1„ 



En 



(16) 



Note that degenerated states are explicitly excluded 
from the sum and therefore |(™|j z (g)|7i)| = and the 
conductivity reduces to a = —n{K)S(uj) that saturates 
the f-sum rule. Interactions spoil the commutation of 
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spin-current and Hamiltonian but, if the Umklapp part 
of the interaction is irrelevant and the system remain 
gapless, the Drude weight remains finite and the current- 
current correlation functions can reduce the Drude peak 
from the kinetic energy. This situation is indeed realized 
in the gapless regime of the XXZ chain but the integrable 
nature of the interaction in this case plays, as we will see, 
a definitive role. The regular part of the conductivity is in 
any case ( integrable or non-integrable systems) enhanced 
to fulfill the f-sum rule. 



a.) 



- s x . ; 













b.) c.) 
FIG. 1. a.) Two-loop contribution to the current-current 

correlation function, a S + and S~ operator must be applied 
in each loop to close it consistently in terms of the loop ori- 
entation, b.) and c.) Two kinds of one loop contributions to 
the current-current correlation function The ordering of the 
S + and S~ in terms of loop time is crucial to evaluate the 
contribution of these terms. 



III. RELATION BETWEEN CORRELATION 
FUNCTIONS ATT/0 



S(q,OJ n ) = —^{[[H, S*], S* ] ) 



4sin 2 ((7/2) 



A(<?,w„) , 
(19) 



where we have used the definition (^). The double com- 
mutator in the right hand side of Eq. ([H]) is the bound- 
ary term of the partial integration and is evaluated to 

([[H,S z q },S*_ q ]) = Asm 2 (q/2)(K) (20) 

Recalling the definition of D(q,w n ) we arrive to 



D(q,u n ) = 



4suT(g/2) 



(21) 



Note, that the double commutator in (|19|) occurs for the 
Matsubara correlation functions and does not occur for 
a related real- frequency correlation functions . 



L=512, T=0.004J 




0.2 0.4 0.6 0.8 

FIG. 2. D(q,u> n ) as a function of u n for various momenta 
q, T — Q.QQAJ XX and J z /J xx = 0.5, 1.5. The lines are fits 
by Eq. (^). Statistical errors are of the order of the symbol 
size. Note the different limiting behavior of D(q — > 0, u>„) for 
Jz/Jxx = 0.5 and J z /Jxx = 1.5. 



Now we will derive a connection in between the 
current-current correlation function D(q, uj n ) and the dy- 
namical susceptibility S(q, u> n ); In the next Section wc 
will explain how to exploit this connection in QMC cal- 
culations. The spin-spin correlation function is defined 
by 

S(q,u n ) = ij\^(S z q (r)S z _ q (0)). (17) 

In Fourier space, the continuity equation (^|) takes the 
form 

^S*{t) = [H,S* g ] = i(l-e*)jf . (18) 

We integrate the right-hand side of Eq. (|l7]) with respect 
to t twice, use (IL8|) and obtain 



IV. QMC EVALUATION OF THE 
CONDUCTIVITY 

In this Section we will discuss the usefulness of Eq. 
( pl| ) in the context of QMC-simulation, comparing two 
different possibilities to compute the conductivity using 
quantum cluster algorithms. 

Cluster algorithms for QMC-simulations allow for 
global updates of the configuration by flipping simulta- 
neously spin-clusters whose typical sizes are of the or- 
der of the_jcorrelation length of the system. The loop 
algorithmEll we used in the present study gives an efficient 
prescription to construct clusters. The resulting autocor- 
relation time is in general of the order of one Monte Carlo 
step (seec3 for an excellent review) . 

The current-current correlation function in real space 
and imaginary time takes the form 
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A(l,r) = 



LNr, 



V,t' 



(22) 



where is the number of Trotter slices. The con- 
tributions to A(l, t) are non-diagonal four-site oper- 
ators, typically (J xx /A)S+(t 1 )S~ +1 {t 1 )S 12 (t 2 )S+ +1 (t 2 ). 
In principle non-diagonal operators can be computed us- 
ing the loop algorithmcJ. When these non-diagonal op- 
erators are two-point-like only one-loop terms contribute 
to the correlation function. In that case it is possible 
to design efficient improved estimators, meaning that a 
given magnitude is evaluated not only in one configura- 
tion but in all configurations related by loop flippings. 
The evaluation, of a four-point correlation function is 
more involvecHj. In that case there are two-loop terms 
and one-loop terms which contribute in different ways de- 
pending on the specific shape of the loop, see Fig. [I] for 
an illustration. As a consequence the improved estima- 
tors are much less efficient. The dynamical susceptibility 
in S(q,u> n ) is, on the other hand, a two-point diagonal 
operator that can be evaluated efficiently using improved 
estimators and it is related to the conductivity using the 
Kubo formula ( |l0| ) and the relation (pl|). 

In particular one can compute within each loop a the 
magnitude: 



W(q,uj n ,a) = 



JL 

N T 



i(ql-\-uj n T) 



(23) 



800 - 2000 and on the average 6 x 10 6 full MC-updates in 
the grand-canonical ensemble. Test-runs within a canon- 
ical ensemble were also performed to exclude any influ- 
ence of the ensemble in the transport properties. The 
error bars (either the statistical ones or those derived 
from the fitting) are of the order of the symbol size in all 
the figures presented. 




FIG. 3. Di(q), Ai(q) and 71(g) from Eq. ( [27] ) as a function 
of momenta q for the XXZ-model, L = 512 and for various J z 
at T = 0.004 J xx . 71(g) is too small for J z < J xx to show up 
on this scale. The lines are the Bethe-Ansatz result (^) for 
the velocity c(J z ) (no fit, for J z < J xx ). For the discussion of 
the fit for J z — 1.5J XX see the text. 



The dynamical structure factor is then 



S(q,w r , 



W{q,to n ,a)W(-q,-Lo n ,a) , (24) 



where a runs over all loops constructed. In particular 
we want to emphasize the importance of relation (21), 



because only using it we obtained the high-quality data 
(large set of uncorrelated measurements with small sta- 
tistical error bars) that is necessary in order to extract 
DC-transport coefficients. 

A second technical important issue is the relation be- 
tween the conductivities in real and imaginary axis. As 
has been discussed by Kirchner, Evertz and Hankc the 
limit lo — > of the conductivity caa be also approached 
from the imaginary frequency axista. Taking the analyt- 
ical continuation to the real frequencies ilo u — > lo + iS 
in the spectral representation of A(q,LO n ) (Eq. |l6|) we 
note that A(q, to + iS) is analytic in the upper half of 
the complex w-plane. Zero-frequency properties like the 
Drude weight or the diffusion constant can be reliably 
extracted by the extrapolation along the imaginary axis 
at low temperatures, when many Matsubara frequencies 
lo„ = 2TtTn are available close to lo = for the extrapola- 
tion. This is the case, however, only at low temperatures 
and will therefore present results only for T <§; J xx . 

Finally we mention a few numerical details. We used 
the discrete imaginary-time version of the Loop Algo- 
rithm with a Trotter decomposition of typically Nt — 



V. DATA-ANALYSIS 

At low temperatures and frequencies, the scaling of 
D(q,uj n ) can be obtained simply invoking the conformal 
symmetry of the model emerging in the gapless regime 
Jz < Jxx- S(q,uj n ) at small q takes then form 



S(q,LO r , 



Di{T)g 2 
(cq) 2 + uj 



(25) 



Note that, unlike near q = ir, the dynamical susceptibility 
about q = do not show power laws. The XXZ-modcl 
maps to an interacting ID spinless fermionic system at 
half filling. For the noninteracting case (the XX chain) 
we can compute exactly D\ [q, u> n ) and we obtain D\ (0) = 
J xx /n, and linv_^ Jxx sin(g) = J xx q = cq. 

Expression J25| ) and Eq. ( pl"| ) suggest the form 



A*(q) 



(26) 



Alternatively, Eq. (26) can be viewed as the first term 



of the exact representation for D(q,u n ) containing an 
infinite-number of termsES: 



D(q,u) n ) 



2 

£ 

3=1 



DM "I 



Af(g) + 27j(g)w„ + w2 



(27) 
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The choice of this fitting function is essential to distin- 
guish the transport properties of ballistic and diffusive 
systems, indeed it allows the correct computation of the 
Drude weight and the diffusion coefficient. We discuss 
now in detail the properties of (p^): 

i. ) D(q,u> n ) is analytic in the upper complex-plane for 
7j(?)>0. 

ii. ) For the zero-q gaps Aj(0) = Um ? _>o A,^) we 
find to possibilities: (a) Ai(0) = and A 2 (0) > 0, 
i.e. Eq. describes a gapless phase, (b) Ai(0) > 
and A 2 (0) > Ai(0) and i.e. Eq. (27) describes a gaped 
phase. In the first case, ( g7j ) reproduces the correct ui 
and q dependence for the scaling form of a the Luttinger 
liquid (pq ) . The first term in Eq. (^) dominates the low- 
frequency behavior in both cases and we have set gener- 
ally 72 = in order to keep the number of parameters to 
a minimum. 

iii. ) At high frequencies 

lim D(0,uj n ) = -(K) = D^+DiiO) (28) 



and a finite D2(0) results in a reduction of the Drude 
weight D{T) with respect to the kinetic energy, see Eq. 
(|12|). A finite D 2 (0) measures therefore the amount of 
decay experienced by the total current due to the in- 
teractions. We note also that the Ansatz Eq. ( f27|) for 
D(q,uj n ), together with Eq. (p8f), is consistent with the 
f-sum-rulec3 



•) 



The uniform spin stiffness p s 



lim^o lim Wji ^ D(q, uj n ) is always zero, as expected for 
a quantum-critical antiferromagnetic chain, 
vi.) The quality of the fit to D(q, u n ) by Eq. ( p7|) is, in 
general, excellent, as illustrated in Fig. ^ for J z — 0.5J XX 
and J z = 1.5 J xx . Note that the q — ► limiting curve for 
D(q,Lu n ) is singular in the gapless phase (J z = 0.hJ xx ), 
but well defined in the gaped phase {{J z — 1.5 J xx ). 
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0.4 



Re <j(ui) duj 



(K) 



(29) 



FIG. 4. The q-dependent Drude weight Di(q) for the 
isotropic Heisenberg chain at T = 0.001J XX for various system 
sizes L = 64,... 512. The T = 0, q = result given by Eq. 
( pTj[ ) is indicated by the arrow. The convergence with system 
size is slow for q — ► 0, due to the logarithmic corrections 
present at the isotropic point J z — J xx . 



for the optical conductivity. 

iv.) In the gapless regime D(q, ui n ) can describe a nor- 
mal conductor with finite DC-conductivity. The optical 
conductivity ( |l0| ) takes for small frequencies the Drude 
form. 



Recr(w) 



2£>i(0)7i(0) 



^ 2 +4 7l 2 (0) ' 1 + (lut) 2 
where we introduced the DC-conductivity 
ao = Di(0)/(2 7l (0)) 
and the quasi-particlc lifetime 

r = (2 7l (0))- 1 . 



(30) 



(32) 



For r — > oo Eq. (]3(]|) reduces to Recr(cj) = nDi(0) 6{w). 
Even more, if we consider now the q-dependence for 
1/ui 3> r, the optical conductivity takes (for small 
c< ?/7i(fJ)) the diffusion form 



a(q,uj) 



a uj 



iD s g 2 



27i (0) 



c 2 r . 



(33) 



D s is the spin-diffusion constant. Eq. ( |33| ) is consistent 
with D s = cX s , where A s = cr is the mean free length. 



VI. BALLISTIC TRANSPORT 

In this Section we will apply the procedure described 
in Section IV to the XXZ chain. We will compare with 
exact known results and study the controversial finite 
temperature behavior of a(oj = 0) for this model. 



In Fig. ||we show the values for D\{q) and Ai(g) and 
7 i(g) for J z / J xx — 0.5, 1.0 and 1.5. These are the values 
used in Fig. |2| for fitting D(q,uj n ). We have included (no 
fit) in Fig. pTfor the gapless regime J z < J xx , the Bethe- 
(31) Ansatz result lim^o Ai(g r ) = c( J z ) q for the magnon dis- 



persion, where c(J 2 ) is the velocity 



ir 8^(6*) 

2 e 



(34) 



2 arcos(J z /J xx ) 

of the des Cloiseaux-Pearson spectrumil, with J z = 
cos(9)J xx . In the gaped phase we have fitted A^g) by 
e{q) = v/A^ + (cq) 2 . We find A = O.igiJ^^j-^vhich is 
close to twice the one-magnon gap of 0.091 J^cll 

The damping 7 i (q) is vanishing small for J z < J xx and 
acquires a finite value in the gaped phase. We found phc- 
nomenological that "fi(q)Ai(q) w const., independent of 
q, for J z > J xx . In Fig. [| we present values obtained by 
QMC for the g-dependent Drude weight for J z = J xx at 
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T = 0.004^2;. We find good convergence for small but 
finite-q, but slow convergence for q — > as a function of 
system size, due to the multiplicative logarithmic correc- 
tions present at the isotropic point. We have indicated 
by thp-horizontal arrow the T = 0, q = Bethe-Ansatz 
resultO. For J z < J xx the agreement in between low- 
T QMCpand the T = 0, q = Bethe-Ansatz result is 
excellentE£l. 

We study now the behavior of the Drude weight at fi- 
nite temperatures for models free from strong multiplica- 
tive corrections. The-jnain conclusion of a Bethe-Ansatz 
calculation by ZotosEa is a fast decay of the Drude weight 
when the temperature increases, in agreement with ex- 
act diagojualization studies in the limit of, infinite tem- 
peratures! Kliimper et al. have foundEj, with an al- 
ternative Bethe-ansatz approach, a functionally different 
behavior for D(T), see Fig. ||. For a numerical probe of 
D(T) we focus on J z = J xx cos(7r/6) and consider several 
small temperatures. For this value of 3 z j J xx the numer- 
ical problems due to multiplicative logarithmic correc- 
tions are absent (compare Fig. ||) and the difference in 
between the two different Bethe-Ansatz predictions are 
substantial. In Fig. (|^) we show a comparipm-pf our 
data with the two available analytical resultsc3£j. Our 
results agree with the temperature-dependence predicted 
by Kliimper et al.. 

We have evaluated the uniform susceptibility x(T) for 
J z = 0.85J XX , L = 512, 1024 and several low temper- 
atures in order to address the two questions: (a) Is it 
correct to compare T = Bethe-Ansatz results with 
QMC-results obtained for a temperature T = 0.004J xa; 
and L = 512? (b) Is T — 0.02J XX large enough for 
L = 512, 1024 not to be affected substantially by finite- 
size effects? The data presented in Fig. ^ shows that 
T = 0.004J.re is indeed below the finite-size gap and 
should be a good approximation to the T — data 
and that for T > Q.012J XX no finite-lattice effect can 
be observed within the statistical error-bars given.. Note 
that the low-T dependence of the Bethe-AnsatzE!l result 
shown in Fig. || for x(T) can be fitted by x(T) ~ T x with 
x ks 0.867. This exponent is very close to the exact value 
x = 0.858 obtained by Eggert et aUB for J z = 0.%5J XX . 
Since x < 1 the slope is diverging for T — > 0. This 

divergence is, however, not relevant for the temperature- 
scale presented in Fig. g. 



0.35 

D(T) 

0.3 
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FIG. 5. QMC results for the Drude weight for 

L = 512, 1024 and J z — J xx cos(-7r/6) = 0.866 Jxx as a function 
of temperature (in units of Jxx)ri*t comparison with two (solid 
lines: Ref.Ej, dashed lines: Ref.Ej) Bethe-Ansatz results. 



VII. PHENOMENOLOGICAL RELATIONS IN 
BETWEEN TRANSPORT COEFFICIENTS 

Let us consider a ID system system with a finite mag- 
netization relaxation time r, which might be due to ei- 
ther intrinsic relaxation processes or due to weak residual 
coupling to an external bath (i.e. phonons). The DC- 
magnetic conductivity takes then the form 



cr(w = 0) = Di(0)T 



(35) 



see Eq. (J3X|) and (]32j). 

Let us assume that an inhomogeneous magnetic field 
B z (x) is applied in the chain along the z-axis. In prin- 
ciple this is a non-equilibrium situation but if dB z /dx 
is small we can assume, after coarse-graining the chain, 
that we have a well defined local magnetization M(x) 
and all thermodynamic relations hold locally. Generaliz- 
ing the usual phenomenologyEil for electric transport we 
write the magnetic current as 



j(x) = v(M(x + A) -M{x)) = Xv 



dM 
dx 



(36) 



where A = vt is the mean-free path v is the velocity as- 
sociated to the magnetization current that we will iden- 
tify latter. We can express the magnetization current in 
terms of the gradient of the magnetic field: 



V 2 T 



dM 
dx 



V 2 T 



dM dB z 
dB z dx 



(37) 



Using j{x) — <J^j— we arrive at the expression 



V 2 TX 



(38) 



This relation is analogous to the well-knownEil phe- 
nomenological kinetic formula for the thermal conduc- 
tivity k = cvv 2 t, where cy is the specific heat and k the 
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thermal conductivity, r can be eliminated if we use Eq. 



a_ = Ji(0)r = D^Q) 
X T XT X 



(39) 



This phenomenological equation is independent of the 
value of t, and holds also in the limit r — > oo, when D\ (0) 
becomes the Drude weight. Since the derivation of ( |39| ) is 
based on quasi-ballistic arguments, it is of interest to ex- 
amine whether this relation holds at low-temperatures 
for Luttingcr liquids and Bethc-solvable models. At 
T = Opbjsth magnitudes x and D have been computed 
exactlyBH for the XXZ chain by Bethe-Ansatz: 



where the first term corresponds to the bosonic part 
and N and J are integer quantum numbers associated to 
states with nonzero charge and current respectively. The 
three velocities present in the Hamiltonian, the sound 
velocity v s , the charge velocity vn and the current veloc- 
ity vj are not independent but restricted by an universal 
relation valid in all microscopic models in the Luttinger- 
liquid universality classE3: 



VNVJ 



(46) 



J xx tt sin(0) 



D(0) = 

4: U I VI — U j 

were we have defined J z — cos(6)J xx and 

x(o) 



Jxxk(it - 0) sin(6>) 



dividing both relations we obtain: 



£>i(0) _ n'J^sm^e) _ 



x(o) 



= c\J z ) . 



(40) 



(41) 



(42) 



For the XXZ chain the values of these three parameters 
of the effective low-energy Hamiltonian ( f45| ) have been 
identified independently, ujv and v s by HaldanecJ and vj 
by Gomez-Santos^ -using the results of different Bethe- 
Ansatz studiesH'E30. 

Guided by the phenomenological derivation presented 
above we propose the following (phenomenological) 
finite-temperature extensions of the velocities in Eq. (|4£ 



1 



v N 



X(TY 



n 



Di(T) 



Jh 



(T) 



c v (T) 



(47) 



This results then allows us to identify the magnetization- 
transport velocity v in (139) with the spin-wave velocity 
c{J z y_ At T = Eq. is then exact. The validity 
of (E2) in leading, low-T correction is an open question 
presently. The leading T-correction to D\(T) and x{T) 
are ~ T 2 for J z — and do not cancel; Eq. (^2|) is ex- 
act for the XX-model only at T = £L The leading T- 
corrections to the susceptibility showa, however, an ex- 
ponent crossover for J z — 0.5J XX and ( p2| ) might hold in 
leading low-T order for J z > 0.5J XX . 

A relation similar |n M (^9[) has been discussed recently 
for jthecmal transportOOexperimentsQ'E3, where we de- 
fineH>EU(T) 



The extension to finite temperatures of vn and vj are 
in agreement with their physical meaning — {^jj)b=o 

and§§@ vj/tt = L(^3^)$ = o for Luttinger liquids. On 
the other hand, for the magnitudes involved in the ther- 
mal ratio, the bosonic part of the Hamiltonian does play 
an important role. In fact, only the bosonic degrees of 
freedom transport the. energy in the homogeneous states 
(which, by dcfinitionEJ do not carry particle currents) rel- 



evant for the thermal conductivity n 



Jh 



t and the spe- 



cific heat cy. It is therefore justified to consider n th /cy 
as the natural extension of v 2 to low temperatures. Re- 
cently this ratio has been computed using Bethe-Ansatz 
techniques at all temperatures by Kluemper and SakaiE-X 
K th j cy is a well behaved function of T, even more it is 
very flat at low temperatures and takes the value the 



(T) t, where K th is-the thermal Drude expected value v s = c (J z ) at T = 0. 



weight. For the XXZ chain one find £3 that both n th (T) 
and cy(T) are linear in temperature for small tempera- 
tures and that 



Jh 



(T) 



lim 

T^O Cy(T) 



Combining Eq. (EM) and Eq. (j43j) we obtain 



£>i(Q) 
x(o) 



Jh 



(T) 



= lim 

T^O Cy(T) 



(43) 



(44) 



For ballistic systems the quantity -Di(O) in above equa- 
tion is identical to the Drude weight. This relation can 
therefore be interpreted in the framework of a Luttinger 
liquid. The Hamiltonian of a Luttinger liquid can be 
written in the diagonal form: 



0.2 
0.18 
0.16 
0.14 
0.12 
X 0.1 
0.08 
0.06 
0.04 
0.02 





/ / 
/ ' 
/ / 



Bethe-Ansatz 

• QWIC, L=512 

• * QMC, L=1024 

T=0, exact 



0.01 0.02 



0.03 



H = J2 Vs \ q \blb g + ±j(v N N 2 +vjJ 2 ) 



2L 



(45) 
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FIG. 6. QMC results for the uniform suszeptibility x(T) 
for L = 512, 1024 and J z = O.&kJxx together with the 
Bethe-Ansatz result (solid line, Ref£j), as a function of tem- 
perature. The star denotes the T — Bethe-Ansatz result. 
Note the absence of finite-size effects for T > 0.012 in the 
QMC data. 



VIII. DIFFUSIVE TRANSPORT 

We study now the effect of non-integrable interaction 
terms in the conductivity of a ID spin system. To be 
specific, we add a small perturbation to H^ xxz \ which 
breaks the integrability of ni xxz )-_ 



H' 



J 



Sa s. 



z+3 



(48) 



The expression (g) for the spin current remains valid, 
since H' does contain only S^-operators and the system 
remains non-frustrated and free from sign problems. We 
have performed QMC simulations for the resulting model 
H = H (xxz ^ + H' mainly for J z = J xx cos(tt/6). We find 
a transition to a gaped phase around J' z = 0.3 J xx , see 
Fig. ^. The exponential opening of the gap resembles 
a Kosterlitz-Thouless transition very similar to the one 
present in the XXZ chain at the isotropic point and sug- 
gests that H' is not changing the universal properties 
of H( xxz ) , only shifting the transition point and adding 
the ingredient of non-integrability. We find the relax- 
ation time t = l/(27i(0)) = liiriq-vo l/(27i(g)) to be 
finite within the numerical accuracy (due to finite-g and 
uj n resolution), leading to a finite DC-conductivity in the 
gapless phase. 

We have examined the temperature-dependence of the 
resulting DC-conductivity. Due to our restriction to 
T <C J xx , resulting from the finite-w„ resolution on the 
imaginary axis (see Section IV) we could not examine a 
large enough T-range in order to determine the full In- 
dependence of cr(0). We found for J' z = 0.3 J xx : a(T = 
0.004J ra ) = 13.6±0.9,cr(T = 0.008 J xa; ) = 12.1 ±1.0 and 
a(T = 0.012 J xx ) = 10.1 ±0.8. 

In agreement with our expectation of a diverging DC- 
conductivity in a translational-invariant system we find 
<r(0.008) > cr(0.012). The increase from ct(0.008) to 
<t(0.004) is, on the other hand, only modest, presumably 
due to the finite-size resolution limitation illustrated in 
Fig. §. 




FIG. 7. For L = 512 and T = 0.004J M the QMC-results 

for the gap Ai(0), the relaxation rate 71(0), the effective 

velocity c and the parameter -Di(O) as a function of J' z for 
H (xxz) + H i with Jz = Jxx cog^/e). 



IX. CONCLUSIONS 

We have shown that Quantum-Monte-Carlo simula- 
tions of quantum-spin chains are a powerful tool to ob- 
tain finite and diverging transport coefficients at very 
low temperatures. We have derived an useful relation 
between the dynamical structure factor S(q,u n ) and the 
dynamical conductivity a(q,uj n ), which allows to calcu- 
late 17(9, w n ) to very high accuracy on the imaginary axis. 
For an integrable chain., we support the original sugges- 
tion by Zotos el alnxA of a finite Drude weight at fi- 
nite temperatures and settle a recent dispute regarding 
the functional form of D(T). In addition we present re- 
sults suggesting the absence of ballistic transport (i.e. 
a zero Drude-weight) for a non-integrable model, for 
which we are able to estimate the magnitude of the DC- 
conductivity. We have discussed our result in the frame- 
light of phenomenological relations and Luttinger-liquid 
theory. Connections to recent studies of the diverging 
thermal conductivity of quantum-spin chains was made. 
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